Localization of resonance eigenfunctions on quantum repellers 
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We introduce a new phase space representation for open quantum systems. This is a very powerful 
tool to help advance in the study of the morphology of their eigenstates. We apply it to two different 
versions of a paradigmatic model, the baker map. This allows to show that the long-lived resonances 
are strongly scarred along the shortest periodic orbits that belong to the classical repeller. Moreover, 
the shape of the short-lived eigenstates is also analyzed. Finally, we apply an antiunitary symmetry 
measure to the resonances that permits to quantify their localization on the repeller. 
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Open quantum systems with classical chaotic behavior 
have been far less explored than their closed counter- 
parts, and as a consequence several of their properties 
are less known. This situation is in contradiction with 
their importance in many areas of physics, some of them 
of very recent development. One of the most significant 
examples is the study of quantum to classical correspon- 
dence ll[ . But their applicability range extends to micro- 
lasers [2, S 3 ' quantum dots [&] , chaotic scattering [1] , 
and more topics. The lack of unitarity at the quantum 
level, and the complex - and mostly fractal - nature of 
the repeller sets is at the heart of the difficulties in their 
understanding. 

The classical phase space of open chaotic systems is 
characterized by the forward and backward trapped sets, 
i.e., fractal sets of trajectories that stay trapped forever in 
the future and in the past, respectively. The intersection 
of both sets is called the repeller. Classical properties of 
this set are its decay rate 7c/, various generalized dimen- 
sions and the spectral properties of the Perron- Frobenius 
operator On the other hand, the quantum evolution 
is described by a nonunitary operator, whose right and 
left eigenstates (the so called resonances) are nonorthog- 
onal and the eigenvalues \i are complex with modulus 
less than or equal to one {uf = \Xi\^ = exp(— F^) < 1, 
where Fi > is the decay rate). The short-lived res- 
onances (Fi ^ 1) are associated with the trajectories 
that escape from the system before the Ehrenfest time, 
while the long-lived ones (F^ = 0(1)) are related to 
the classical forward or backward trapped sets (in 
they were coined quantum fractal eigenstates). Accord- 
ing to the conjectured fractal Weyl law, the number of 
long-lived resonances can be related to h and the 
structure of the classical phase space through the ex- 
pression N-y ^ [d is the fractal dimension of the 
classical strange repeller). This law has been checked 
for a three disk system Q and some quantum maps 

i, [13, [HE E Qui. 

Finer information about the semiclassical limit is pro- 
vided by the structure of eigenfunctions, about which 



much less is known. For the open baker map it was found 
that the probability density averaged for several right 
eigenstates is supported by the corresponding classical 
forward trapped set p^] . A simple argument by analogy 
to Shnirelman's theorem would favor a weak equidistribu- 
tion of the single spectral components over the repeller, 
but of course this leaves ample room for "accidental" lo- 
calization on the periodic orbits, i.e. for the existence of 
scars in open quantum systems. Such long-lived scarred 
modes have been observed in optical microcavities [4|, 
and have been extensively studied in closed systems [16 1 . 
Recently, some results for the open cat map have been 
provided in [14]. But the knowledge we have about this 
behavior is still rudimentary. 

In this letter we show that many long lived resonances 
of the open baker map, besides the obvious localization 
on the repeller, indeed show an enhancement on some of 
its shortest periodic orbits. To this purpose we introduce 
a phase space representation of the spectral components 
of the open map that is in complete analogy to the usual 
Husimi functions of the closed case. We also investigate 
the shape of short-lived resonances of the open triadic 
baker map. We find that structures resembling the states 
of a simple Hamiltonian develop outside of the classical 
repeller. On the other hand, we have applied a measure 
of antiunitary symmetry to the eigenstates. By doing 
this we were able to quantify the localization of a given 
resonance on the repeller. Moreover, this measure turns 
out to be simply related to the moduli of the eigenvalues. 
We believe that these tools could provide an alternative 
way of identifying long and short-lived resonances, a task 
of the utmost importance in the road to demonstrate the 
conjectured Weyl law In fact, this stands as an open 
problem, where only for the Walsh quantized open baker 
map [H, [l3| exact results could be obtained. 

We consider a closed quantum map U, chosen as the 
quantization of a classical map which is chaotic. The 
phase space is the unit torus and upon quantization we 
impose antiperiodic boundary conditions resulting in a 
Hilbert space of dimension N, where h = {2ttN)-^ plays 
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the role of an effective Planck's constant. U is then an 
N X N unitary matrix. Our maps have parity (i?) and 
time reversal (T) symmetries in both its classical and 
quantum versions. 

The map is "opened" as usual by means of a projection 
operator Ho that defines the escape region in phase space. 
In this way, 11 = 1 — Hg composed with the closed map 
defines the open system. Depending on the symmetries 
of the projector it is possible to maintain or break in 
the open map the symmetries of the closed one. The 
resulting operator, U = UH is not normal and has a null 
subspace whose dimension is tr{Ilo). If we assume that 
it is diagonalizable, its spectral decomposition in terms 
of right and left eigenstates is 



U 



(1) 




Figure 1: Autocorrelation function in phase space for the 4th 
iteration of the standard quantum baker map for = 160. 
The closed version, showing all periodic points of period 4 
is on the left panel, and the open version with i = 3 is on 
the right panel. The periodic orbits which have the forbidden 
sequence (000 and 111) are absent in the open system. 



where we assume the normalization condition {ipiltpi) = 

= 1. As the map is contracting, the spectrum 
is strictly inside the unit circle. Using the decomposition 
in Eq[T]the autocorrelation of a coherent state \q,p) can 
be written as 



{q,p\U''\q,p) 



in terms of complex phase space distributions 



hi{q,p) 



{q,p\^^){lP^\q,p) 



(2) 



(3) 



As the right (left) eigenfunctions are approximately con- 
centrated on the backwards (forward) trapped sets, we 
expect \hi{p,q)\ to map out the distribution of the sin- 
gle resonances on the intersection of these sets, i.e. 
on the repeller. Thus they constitute the ideal tool 
to study the localization properties and eventually the 
"scarring" in open systems. Notice that \hi{p,q)\ (x 
^/nfip.q) n'^ip.q), where Uf ' are Husimi functions of 
the i?, L eigenstates, and that for unitary systems they 
reduce to the Husimi functions of the eigenstates. We 
would like to mention that, instead of considering coher- 
ent states we could have chosen any other representation 
in phase space. Different choices can help to analyze dif- 
ferent properties of \ipf'){'fpi\/ {'>pf\i^f') ■ 

We apply these ideas to study the localization of 
resonance eigenfunctions for the open baker map. In 
its standard version it is defined as {p' = p/2 + 
\2q\/2,q' — 2q — [2gJ), where [2q\ is the integer 
part. It has a complete symbolic dynamics in terms 
of two symbols. The quantization is obtained in 
terms of the anti-periodic, discrete Fourier transform 
(Gjv)j-,fe = expi-i2TT/N(j + l/2){k+l/2)/^/N, as U ^ 

G]y diag{GN/2,GN/2) [3, E^. The advantage for this 
kind of study is that one can easily engineer openings 
that coincide with some iteration of the basic Markov 



partition and thus obtain repellers with a simple sym- 
bolic description in terms of a subshift of finite type. 
For concreteness we take the escape region Hq to be the 
union of the areas defined by qi G [0, 1/2'] x pi e [0, 1] 
and 92 G [1 — 1/2', I] x p2 & [0, 1] leading to a subshift 
for which the complete symbolic dynamics is pruned of 
all trajectories that include the strings O' and l'. The 
opening is symmetric with respect to q = 1/2 in order 
to retain the parity symmetry. We will focus on the first 
non-trivial case oi I = 3 where qi E [0, 1/8] q2 G [7/8, 1], 
and the strings 000 and 111 are pruned. The topological 
entropy At ~ 0.4812 is easily calculated from the leading 
eigenvalue of the transition matrix. In Fig. [T] we com- 
pare the autocorrelation function in phase space (Eq. [2|) 
for both the open and the closed versions of the fourth 
iteration of the map clearly showing the disappearance 
of the periodic points forbidden by the pruning rules. 

Four examples of resonance eigenfunctions are dis- 
played in Figl21 using the distribution hi{q,p). They 
correspond to a Hilbert space of dimension = 320 
and 1 — 3. In panel a) we show the longest-lived reso- 
nance for which |Aa| ~ 0.977. In this case, the complex 
function h{q,p) is represented through an intensity (mod- 
uli) and color (phase) scale, this latter ranging from blue 
(0) to red {2tt). It is clear that its values are almost 
completely real, this being an expected property of an 
eigenstate whose associated eigenvalue is so close to the 
unit circle. In panels b), c) and d) we show \hi{q,p)\ for 
eigenstates corresponding to |A;,| ~ 0.905, |Ac| « 0.856, 
and \Xd) \ ~ 0.641 respectively (in these cases we restrict 
our analysis to the moduli of the distribution) . We now 
focus on the specific shape of resonances. Those in a) 
and c) clearly show an enhanced probability on the peri- 
odic orbits 0011 and 01 respectively, which belong to the 
repeller. This enhancement on the periodic points and 
their manifolds is typical of the scarring phenomenon. In 
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Figure 2: (color online) Phase space representation of a set 
of resonances for the standard open baker map with A'' = 320 
and Z = 3. In panel a) we display hi{q,p) for the longest-lived 
resonance (|Aa| « 0.977). Complex phase is represented by 
color ranging from blue (0) to red {2n). In b), c) and d) we 
show \hi{q,p)\ for eigenstates corresponding to |A(,| « 0.905, 
|Ac| ~ 0.856, and \Xd)\ ~ 0.641, respectively. The classical 
repeller for times between t = ~2 and t = 5 lies inside the 
solid red line rectangles. 



contrast, b) has no visible prevailing orbit in its structure, 
and it is distributed almost uniformly on the repeller. We 
have represented a finite time version of the classical re- 
peller by means of red line rectangles as a guide to the 
eye. We have taken times between t — —2 and t = 5 to 
keep symmetries intact. Finally, in order to complete the 
picture of a representative set of resonances, we include 
one that it is localized out of the repeller, in panel d). 

In order to study the eigenstates outside of the repeller 
in a clear way, we choose to analyze the triadic baker 
map Bt- This map is defined as {p' = p/3,q' — 3q) if 
< q < 1/3, (p' = (p+l)/3, q' = 3g-l) if 1/3 < Q < 2/3, 
and (j)' ^ {p + 2)/3, q' = 3q - 2) if 2/3 < g < 1. From 
our point of view, the main advantage of this map is that 
it has an isolated fixed point at {q,p) = (1/2,1/2), far 
from the discontinuities. The quantization of this map in 
a TV-dimensional Hilbert space {N = 27r/fi,, divisible by 
3) with antisymmetric boundary conditions is given by 
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Figure 3: (color online) Phase space representation \hi{q,p)\ 
of the six longest-lived resonances for the triadic open baker 
map with A'^ — 243. Moduli of the corresponding eigenvalues 
are \Xa\ = l/VS ^ 0.577, lA^I = |Ac| « 0.249, |Ad| = |Ae| ^ 
0.184 and |A/| « 0.064. The map is real, therefore \a,Xf € R, 
Af, = A*, and Ad = A^. 



The open version of the map can be obtained in the same 
way as before, by means of a projector 11, but this time 
on the central region in position. This amounts to saying 
that in the quantum case the left and right thirds of the 
columns of Bt are set to zero. This opens two thirds of 
the phase space (where the forbidden region is given by 
q e [0,1/3] and q G [2/3,1], satisfying tr(Il) = N/3 in 
the quantum case) . The repeller of this system is reduced 
to the point {q,p) = (1/2,1/2). It is important to un- 
derline that this map is hyperbolic, but it is not chaotic. 
However, it is very suitable in order to unveil the struc- 
ture of resonances lying outside of the repeller, that in 
the previously studied version of the baker map have a 
complicated structure. Fig. [3] shows the six longest-lived 
resonances for iV = 3^. In this case the morphology is 
very clear. They resemble excited states of the hyper- 
bolic Hamiltonian H = qp, centered at the fixed point 
{q,p) = (1/2,1/2). This suggests that the longest-lived 
fraction of the short-lived resonances associated to more 
complicated repellers could be interpreted (and perhaps 
theoretically constructed) as convenient superpositions 
of these simpler structures. This completes a very inter- 
esting picture, i.e., one could try to devise a theory of 
short periodic orbits capable of describing the complete 
significant quantum information of an open system (20j . 

Finally, the time reversal symmetry provides us with 
very interesting results regarding the classification of res- 
onances into long and short-lived. In the classical case 
this symmetry can be restricted to the orbits belonging to 
the repeller, since they are not affected by the presence 
of the escape regions, in our example. In other words, 
the repeller corresponding to our open system retains 
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Figure 4: (color online) Measure of time reversibility of a given 
eigenstate (n = \{ipf'\T\'ipi)\) as a function of the modulus of 
its eigenvalue, for the standard open baker map with / = 
2, 3, 4, 6 and A'^ — 1024. Dashed (gray) lines correspond to 



spending classical repeller. Moreover, we have obtained 
some hints on the structure of the short-lived eigenstates, 
i.e. those living outside of the invariant set. We believe 
that short periodic orbits of the repeller play an impor- 
tant role also in this case. Finally, a time reversibility 
measure for the resonances has been devised. We were 
able to find a simple relation with the eigenvalues moduli 
showing that the long-lived resonances have the highest 
time reversibility. We consider that these tools are very 
promising in the way to demonstrate the fractal Weyl law 
for open quantum systems. 

Partial support by CONICET and ANPCyT is grate- 
fully acknowledged. 



the time reversal invariance of the closed one. However, 
if we look at resonances in Fig [2] we immediately notice 
that in general they are not symmetric with respect to 
the line q = p, this indicating that they are not time 
reversal invariant (in general left eigenstates cannot be 
obtained from applying the time reversal operator to the 
right ones, for example). Then, for the quantum case 
we can define a measure of the time reversal symmetry 
for each eigenstate in order to see if this is related to 
its localization on the repeller. This is accomplished by 
n EE (V-f iTlV'f ), where T = KG {K stands for the com- 
plex conjugation operator). In the case of maps having 
a time reversal symmetry ti — 1 for any i. Results for 
the standard open baker map can be seen in Fig. [4l We 
have found that there is a simple relation among the val- 
ues of this measure and the eigenvalues moduli, given by 
Ti = |A|'~^. On one hand this confirms that the reso- 
nances that have a maximum overlap with the repeller 
(long-lived ones) show the higher time reversibility. On 
the other hand, it confirms this quantity as an alterna- 
tive way to distinguish between long and short-lived res- 
onances. We think that by following this approach a pre- 
diction on the number of each kind of eigenstates could 
be found 

In conclusion, we have introduced a very powerful 
phase space representation for open quantum systems 
that is in complete analogy with the usual Husimi distri- 
bution for closed ones. This allowed us to clearly show 
that long-lived resonances can be strongly scarred along 
the shortest periodic orbits that belong to the corre- 
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